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Vortex Dynamics in Ferromagnetic Superconductors: Vortex Clusters, Domain Walls and 

Enhanced Viscosity 
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We demonstrate that there is a long-range vortex-vortex attraction in ferromagnetic superconductors due to 
polarization of the magnetic moments. Vortex clusters are then stabilized in the ground state for low vortex 
densities. The motion of vortex clusters driven by the Lorentz force excites magnons. This regime becomes 
unstable at a threshold velocity above which domain walls are generated for slow relaxation of the magnetic 
moments and the vortex configuration becomes modulated. This dynamics of vortices and magnetic moments 
can be probed by transport measurements. 

PACS numbers: 74.25.Uv, 74.25.F-, 74.25.Ha, 74.25.N-, 75.60.Ch 



Introduction - Superconductivity (SC) and magnetism are 
at the heart of modern condensed matter physics. While 
they seem to be antagonist according to the standard BCS 
theory, a large family of magnetic superconductors were 
discovered in the last decades. Examples include coexis- 
tence of antiferromagnetism or helical ferromagnetic (FM) 
order in ternary superconducting compounds! D> uniform fer- 
romagnetism in triplet superconductors! 2-4 1, an d antiferro- 
magnetism in the ReNi2B2C borocarbides[5 1 (Re represents 
a rare earth element) and in the recently discovered iron- 
based superconductors|6|. The interplay between SC and 
magnetism allows to control the superconducting properties 
through the magnetic subsystem, and vice versa. These phe- 
nomena open new possibilities for applications to supercon- 
ducting electronics and magnetic storage devices Q 0. 

The Abrikosov vortices of superconductors are a natural 
link between the superconducting condensate and the mag- 
netic moments. Vortices are induced either by external mag- 
netic fields or by the MMs |9 |. On the other hand, the mag- 
netic subsystem supports collective spin-waves and topolog- 
ical excitations that are domain walls. Because vortices are 
magnetic objects, they are expected to interact strongly with 
MMs via Zeeman coupling. Indeed, as we discuss below, vor- 
tex motion can drive magnetic domain walls. 

The MMs provide a novel handle to control the vortex be- 
havior in the static and dynamic regimes. It was demonstrated 
that magnetic domains induce a vortex pinning that is 100 
times stronger than the one induced by columnar defects lflOl . 
In the flux flow regime, vortex motion radiates magnons by 
transferring energy into the magnetic system. This effect has 
been recently proposed by Shekhter et al. for antiferromag- 
netic superconductors ifTTI . By assuming a rigid vortex lat- 
tice and fast relaxation of the MMs, it is demonstrated that 
Cherenkov radiation of magnons occurs when the vortex lat- 
tice velocity, v, satisfies G • v = £2(G), where G is the vortex 
lattice wave vector and £2(G) is the magnon dispersion. This 
emission gives an additional contribution to the vortex viscos- 
ity that manifests as a voltage drop in the I-V characteristics. 
Thus, the overall dissipation is reduced for a given current. 
Vortex motion can also be used to probe the spectrum of exci- 
tations in the magnetic subsystem. lfT2l 



Several questions remain to be addressed. It is known that 
intrinsic nonlinear effects of the magnetic subsystem become 
important for high energy magnon excitations. However, it is 
unclear if magnon excitations remain stable in this nonlinear 
regime. On the other hand, the interaction between the mag- 
netic subsystem and vortices may become comparable or even 
stronger than the inter-vortex repulsion. Therefore, the vortex 
lattice may be modified by this effect. Finally, the dominant 
dissipation mechanism of vortices when domain walls are ex- 
cited by the vortex motion is unknown. 

Here we study the vortex dynamics in FM superconductors. 
The Zeeman coupling between vortices and MMs induces an 
additional vortex-vortex attraction that is comparable to the 
inter-vortex repulsion for a large enough magnetic suscepti- 
bility. This attraction leads to the formation of vortex clusters 
at low vortex densities. We also show that magnetic domain 
walls are created when vortex clusters driven by the Lorenz 
force reach a threshold velocity. The interaction between do- 
main walls and vortices greatly enhances the vortex viscosity 
and causes hysteresis in the dynamics of the whole system. 
The vortex configuration is modulated by the domain walls. 

Model- Uniform FM order and SC suppress each other be- 
cause of the exchange and electromagnetic coupling between 
the MMs and Cooper pairs [1 J. However they could coexist in 
triplet FM superconductors ||2]-@), such as UGe2, layered mag- 
netic superconductors consisting FM and SC layers lfT3lfT4l . 
such as Smi 85Ceo 15G1O4, or artificial bilayer systems ll8l [T5l . 
Here we study the vortex dynamics in these FM supercon- 
ductors. An applied dc magnetic field perpendicular to the 
ferromagnetic easy axis creates a vortex lattice that is driven 
by a dc in-plane current. We use the approximation of straight 
vortex lines and the description of vortices is two dimensional. 

The total Gibbs free energy functional of the system, in 
terms of the vector potential A, magnetization M and vortex 
position Ri = (x,,y,), is 



G(A, M, RO = d j 



dr (gsc +gM + gint)+ 



1 f 

8tt J ou , 



dr 3 B 2 , (1) 



where d is the thickness of the system and the last term is 
the magnetic energy outside the superconductor. The energy 
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functional density for the SC subsystem in the London ap- 
proximation is 



Ssc(A) = 



B B ■ Hp 



1 



An 



V0-A 



(2) 



with B = V x A. </> is the superconducting phase, H ext is 
the applied magnetic field, Ai is the London penetration depth 
and On = hc/2e is the flux quantum. The energy functional 
density of the magnetic subsystem is 



gM = J ( VMf- J 4M 2 x , 



(3) 



where J and J a are the exchange and anisotropy parameters. 
The easy axis is taken along the x direction. We assume 
that the magnitude of the magnetic moment is conserved, 
\M\ = M„ where M s is the saturated magnetization value. 
Because of the anisotropy, the magnetic Hamiltonian has two 
degenerate minima and supports stable domain walls. The 
Zeeman interaction between MMs and SC is 



= -B M. 



(4) 



The vortex axis is taken along the z direction. The straight 
vortex lines approximation is valid when d <K Ai or d » A^. 
The spreading of magnetic field associated with vortices near 
the surface of superconductors has to be taken into account for 
d ~ Ai, [16|. By minimizing g sc + g mt with respect to A, we 
obtain the magnetic field associated with vortices 

x V x (B - 4/rM) + B = O V 6(r - R;)l (5) 

i 

M z (k) = B z (k)x zz (k) in the linear response region when 
M z /M s <k 1. As Ai is much larger the magnetic correla- 
tion length £ m ~ -\PJ7a, we can use a local approximation 
f° r Xzz(k) — 1/Ja = aW(1 + 4^fo)- The uniform susceptibil- 
ity Xo K (M z (h = 0)M z (k = 0)) diverges at J A = An, which 
signals an instability of the magnetic subsystem. The FM or- 
dering along the x-direction coexists with superconductivity 
only when J A > An. [17] 

According to Eq.|5|, the magnetic field of a vortex at R, is 



B z (k,Ri) = 



On 



1 + Aik 2 



exp(/k ■ R, ), 



(6) 



with a renormalized penetration depth A e = A L / Jl + Anxo- 

Attraction between vortices via MMs - We calculate now the 
interaction between two vortices at R, and R ; . Vortices inter- 
act with each other through the exchange of massive photons 
described by g sc , which leads to a short-range repulsion. As 
was first discussed by Pearl, vortices also interact through the 
exchange of massless photons outside the SC, as described by 
the last term in Eq. (QJ. This contribution leads to a long-range 
repulsion] 18 19]. The total repulsion energy is 



U r (R) 



8n 2 A 2 



K 



8nA 



He, 



(x)-(x) 



(7) 




FIG. 1. (color online) Vortex-vortex interaction potential for different 
values of^o according to Eqs. |7| and {8}. Attraction is induced due 
to the Zeeman coupling between vortices and MMs, and the long- 
range repulsion arises from the electromagnetic fields outside the SC. 



with R = R, ■ - Rj and A = 2A e coth[d/A e ] is the modified 
Pearl length. Kj is the modified Bessel function, Hq is the 
Struve function and To is the Weber function. 

A vortex at R ; polarizes the surrounding MMs. This effect 
leads to an effective attraction to a vortex at Rj. The magnetic 
energy due to the presence of vortices is d J dr 2 (gM+gmt) with 
Bl = B z (Ri) + B z (Rj) and M\ = xoB l z /(l + Anxo). The contri- 
bution from the gradient term in Eq. Q is much smaller than 
the anisotropic contribution because k£ m <K 1 with k ~ 1/A e . 
By using M\ + M| = M 2 S , we obtain the attractive interaction 



U a (R) 



I 



drEtMi 



An(l +4nxo)*> 



\A„ 



(8) 



In the presence of attraction, the repulsion through the elec- 
tromagnetic fields outside the SC in Eq. |7| cannot be ne- 
glected because it prevents the formation of a single cluster. 
The physics here is similar to the laminar phase in conven- 
tional type I superconductors [20 1 . 

The effect of finite velocity, v, on the vortex-vortex in- 
teraction is negligible because Xzz depends weakly on v for 
g m /A e <k 1. The attractive component is comparable to the 
repulsion for xo ~ 1 and the energy minimum takes place at 
R m ~ A e . Fig. ([TJ shows the energy of two vortices separated 
by a distance R. For xo ~ L the net interaction is attrac- 
tive for large separations A e < R < A and repulsive at short 
distances R < A e . There is also a long-range repulsion for 
R > A due to the surface effect. Since the susceptibility xo de- 
creases with J a, the attractive component drops as anisotropy 
increases. The inter- vortex interaction becomes purely repul- 
sive for^-o <K L 

Excitation of domain walls - We introduce the equation of 
motion for MMs and vortices that is used in the numerical 
simulation. The FM subsystem is described by the Landau- 
Lifshitz -Gilbert equation BTII 



d,m = -ym x B e ff + am x d,m, 



(9) 



where y is the gyromagnetic ratio, m = M/M s is the nor- 
malized MM, a is the damping coefficient and the effective 
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(a) / e xt,y = 0, E y = (b) J exty = 0.08, E y = 0.308 (c) / ext , y = 0.1, E y = 0.473 (d) / exty = 0.12,% = 0.604 




FIG. 2. (color online) (a-b): Development of the amplitude of superconducting order parameter |*P|, magnetic structure (x-component of 
the magnetic moment: M x ) and B- as the current increases. The vortex positions correspond to the regions with suppressed (a) Static 
configuration with 7 ext y = 0. In (b-d), domain walls are created and the vortex configuration is modulated. |*P| is suppressed (top row) and B- 
is maximal (bottom row) in the normal core of vortices. MMs are canted by vortices so M x is reduced (middle row). 



magnetic field is B e ff = S[gM + gint]/(5M. The vortex sub- 
system is described by the time-dependent Ginzburg-Landau 
equations 



2mD 



d, x ¥ = - 



a*¥ + B m 1 ¥ + — I iV + —A 
2m \ (D 



, (10) 



become decoupled. Therefore, we shall consider the interest- 
ing region B z < B s . 

In the long wavelength and weak damping a <K 1 limits, 
the magnon dispersion for the FM system of Eq.|9]) is 



Q 2 = u>l + v 2 s k 2 , v s = yM sy l(2-m 2 z0 )J A J, (13) 



-d,A = J s + J ext - -pV x (V x A - 4ttM), 
c An 



with the supercurrent 

eti 



4e 2 



(^•v^ - w^*) - — m z A, 

im mc 



(11) 



(12) 



D is the diffusion coefficient, cr is the conductivity in the nor- 
mal state, Je X t is the external current and other parameters are 
defined according to the usual convention. The MMs stop 
responding to the vortex motion when the average magnetic 
field, B z x n v (£>o with n v being the vortex density, is larger 
than the saturation value, B s = M s Ja, and the two subsystems 



J 2 A y 2 M 2 (l-m 2 z0 )[l+ia(2- 



4r 



r 1/2 ], 



(14) 

where m z o is the z component of the MMs in the ground state 
and v s is the magnon velocity. Re(wo) is the energy gap and 
Im(wo) is the magnon relaxation rate. Re(w ) = 100 GHz and 
v s = 50 m/s for typical ferromagnets. E2l 

We then establish general relations of the energy transfer 
between MMs and vortices. The vortex velocity acquires an 
ac part, v,-, because of the interaction between vortices and 
MMs, v, = v + Vj. The energy balance for the whole system 
reads 

rjv 2 +n(vi 2 ) + -77- { f ^ 2 (d,M) 2 \ = ¥ L ■ v, (15) 
\ n v M s y\J l xl 
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where (• • ■ ) it denotes average over vortices and time, and 
(■ • ■ ) xt denotes average over space and time. The first and 
second term on the left-hand side (lhs) correspond to Bardeen- 
Stephen (BS) damping with coefficient rj = (D ( 2 ) cr/(2^c 2 ^ 2 ) 

GUI , where £ = -Jn 2 /(2m\a s \) is the coherence length. The 
third term on the lhs accounts for the dissipation due to preces- 
sion of MMs. The term on the right-hand side is the work done 
by the Lorentz force Ft,. The effective viscosity r]^ = F^jv is 
enhanced due to the interaction between vortices and MMs, 

/7eff = 7? + \ (vA dr 2 (d,M) 2 ) ■ (16) 

v 2 N 'M n v v z M s y \J J xt 

Off resonance, the contribution of the magnetic damping is 
small, thus v ~ Fl/?]- Since Fi = J ext <&{)/c and E = v/i v Oo/c 
with an external current J ext an d electric field E, the underly- 
ing dynamics can be probed by the I-V measurement. 

The effect of magnons on the vortex dynamics depends on 
the vortex density. When the average inter-vortex distance is 
smaller than the value corresponding to the potential mini- 
mum, n v < l/R 2 t , the attraction between vortices dominates. 
Vortices form circular clusters with internal triangular struc- 
ture in the ground state, as shown in Fig. |2](a) obtained from 
our simulations [23 1 . The distance between neighboring vor- 
tices inside the cluster is of order A e , and the separation be- 
tween neighboring clusters is of order ^nRll{n v A 2 ), with a 
cluster radius given by R c as A[-m /(3m j .)] 1/3 .[24| The attrac- 
tive, u a < 0, and repulsive, u r > 0, energies are defined in 
Fig. [TJ The vortex clusters start to merge and more complex 
vortex configurations, such as stripes, are possible for larger 
values of n v . The H = H c \ transition from the uniform Meiss- 
ner state to the state with vortex clusters is of first order ll25l - 
l27l in contrast to the second order phase transition expected 
for conventional type II superconductors [20]. Vortex clusters 
in conventional superconductors with inter-vortex attraction, 
such as Nb, have been observed experimentally, see Ref. |28| 
for a review. 

For finite transport current, each cluster driven by the 
Lorentz force moves as a whole and polarizes the MMs along 
its way. The MMs relax to their positions of equilibrium 
after the vortex cluster leaves that region. The polariza- 
tion and excitation of magnons, and subsequent relaxation 
of MMs thus causes vortex dissipation through the magnetic 
subsystem! 29 1 . The static structure of the vortex clusters re- 
mains the same for a small v because the change of the vortex- 
vortex interaction is negligible for g m /A e <k 1. 

Here we derive a resonant condition between the motion 
of vortex clusters and magnon emission. The magnetic field 
distribution produced by the vortex motion has a dominant 
wave vector G x = 2n/R m , with R m w A e as shown in Fig. 
[TJ The unperturbed ordered state has M-o = 0. The resonant 
condition G x v = 0(G A ) gives a resonant velocity for vortices 
moving along the x direction 




FIG. 3. (color online) Difference between the electric fields induced 
with and without magnetic moments as a function of current J exUy , 
AE = Em - Eb, where Em is the electric field for the system with 
magnetic moments and E B is the electric field for the system with- 
out magnetic moments. The vortex viscosity increases when domain 
walls are created resulting in a drop of the electric field (vortex ve- 
locity). 



This linear analysis is correct as long as the canted MMs sat- 
isfy the condition that M zc ~ <S> /(J A R 2 n ) •« M s [or J A » 
O /(R 2 m M s )]. 

The oscillation amplitude of MMs and the ac part of the 
vortex velocity are greatly enhanced in resonance and ^ft in- 
creases according to Eq. 



16 1. Two competing processes are 



involved in the magnetic subsystem: the energy input from 
vortex motion and the magnetic relaxation. For large dissi- 
pation (a » 1), the excited magnon is quickly dissipated 
and the vortex cluster with canted MMs remains stable. On 
the contrary, the incoming energy accumulates for weak mag- 
netic dissipation, a «: 1, and increases with time. This 
effect leads to an instability of the magnon excitations that 
has been discussed decades ago both experimentally QUI and 
theoretically 1 3 1 ~ [33ll . For a large enough oscillation ampli- 
tude, the MMs are no longer restricted to one of the symmetry- 
breaking states (there are two degenerate ground states with 

MxQ = ±M S Jl — m^j) and they can flip to the other ground 
state (with opposite M x q). Domain walls are then created as 
shown in Fig. |2|b, c, d). M z becomes large inside the do- 
mains walls and this effect increases the coupling between 
the magnetic subsystem and vortices. For v » v f , the clus- 
ter structure evolves into vortex stripes along the driving di- 
rection Qb, c, d)]. The domain walls are oriented along the 
vortex stripes due to the strong attraction between vortices and 
domain walls. Vortex stripes for large driving forces and ran- 
dom pinning potentials have also been observed in numerical 
simulations without MMs f34l . As vortex clusters drive do- 
main walls, the dissipation increases and the vortex velocity 
(voltage) drops as shown in Fig. [3] The threshold velocity ob- 
tained from simulations where the domain walls are created is 



compatible with that estimated from Eq.( 17 1 



v, = yM s \\2J A J + 



R 1 J 2 
4n 2 ' 



(17) 



Discussions - The magnetic susceptibility is small, %q <sz I, 
in bulk FM superconductors such as UGe? Il35l . Thus, the at- 
traction between vortices is negligible and the ground state 
is a triangular vortex lattice. In the flux flow regime, the 
vortex lattice is resonant with the oscillations of MMs when 
G • v = £2(G) is satisfied. We predict an enhancement of the 
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vortex viscosity at resonance, which can be probed by the I- 
V measurement. A large susceptibility, xu ~ L is needed to 
realize the vortex cluster configuration. This requirement can 
he fulfilled by some cuprate superconductors with rear-earth 
elements (Re), such as ReBa2Cu30 x , where Re ions order an- 
tiferromagnetically below ~ I K. Spins are free from the 
molecular field above the Neel temperature TV ~ 1 K and 
can be easily polarized 061 l37l to mediate the attraction be- 
tween vortices in the low magnetic field region. The vortex 
cluster phase can also be achieved in heterostructures of su- 
perconductors and ferromagnets with large susceptibility 1 38 1. 
On the other hand, random pinning centers may prevent the 
formation of vortex clusters because pinning is strong for a 
small vortex densities. However, vortex motion in the flux 
flow regime quickly averages out the effect of random pinning 
centers 1 39, 40] and the cluster structure may be recovered. 
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FIG. 4. (color online) (a) Approximated inter-vortex potential inside 
a vortex cluster. A s a A e and A, = A compared with the interaction 
profile in Fig. 1 of the main text, (b) Partition of the vortex cluster 
with radius R c to calculate the energy of the cluster. 



Supplement : Characterization of the vortex cluster phase 

Here we characterize the vortex cluster phase by using a 
simple model. The vortex clusters form a triangular lattice 
with lattice constant a c . Each cluster contains n c vortices. In 
the dilute vortex phase, the interaction between vortex clusters 
is Coulomb-like, 



® 2 n 2 c y L_ 



(18) 



where R c j is the position of the cluster and the summation 
is over all clusters. The summation is calculated numerically 
and the result is well approximated by the expression 



An 2 a r 



(l.6AT c 15 -2N C ), 



(19) 



where N c is the number of clusters. N c = n v L 2 /n c and a c = 
Lj y[W c in a sample with lateral size L 2 and vortex density «,.. 
The repulsion between all clusters is 



ttR 2 /A 2 . Since the /^-independent term of Eq.(f20| is irrelevant 
in the following calculations, we will neglect it. 

To estimate the interaction energy inside the cluster we ap- 
proximate the inter-vortex interaction as shown in Fig. |4| a )- 
A vortex in the region with radius R c - A, attracts vortices in 
a region nA 2 and repels vortices in a region n (r 2 - \ fj [see 
Fig. |4|b)]. Then, the interaction energy in the region with 
radius R c - A s is 

Ua = ^ (R c - A,) 2 [nA 2 s u a +n(R 2 c - A 2 .) u r ] , (21) 

where u a < is the attraction and u r > is the repulsion. 
A vortex in the ring R c — A s < r < R c attracts less vortices. 
The attraction region for a vortex in the ring can be written 
as a' A 2 , with or' » 2 obtained by direct integration over the 
ring area. The repulsion region is given by nR 2 — a' A 2 . The 
interaction energy in the ring is then given by 

t/c2 = ^ A, [2A]u a + (ttR 2 . - 2A 2 S ) u r ] . (22) 



The total interaction in the vortex clusters is 



U c = N c (U ei + U a ) ■ 



(23) 



In thick superconductors with d » A, we have \u a \ » u r 
and R c » A s . Then the energy of the whole system (apart 
from the R c -independent contribution) is 



U ~- 



®ln v ^,L 2 R c [2n-K 



-R c A s u a + 2^ R c u r 



n v L z A 



2 )2 



nR 2 ' 
(24) 

where the first term accounts for the long-range repulsion be- 
tween vortex clusters, and the second term accounts for the in- 
teraction inside clusters. The interaction between vortex clus- 
ters scales with the density as nj 2 and the interaction inside 
the cluster scales as n v . The first term can be neglected in the 
dilute vortex case n v A 2 « 1. By minimizing U with respect 
to R c , we obtain the radius of the vortex cluster 



(25) 



Ogn v y/JT v L 2 R c 
2nV~ 2 X 



1.6 



<Eq(^ 2 ) 2 
4tt 2 L 



(20) 



where we have used N c = n v L 2 / (nR 2 /A 2 ). R c is the cluster 
radius and A s is the separation between two nearest vortices 
inside the cluster. The number of vortices in a cluster is n r = 



and the distance between nearest vortex clusters is 

r^T R,. 



a c = .J — — . 



(26) 



By comparing the potential in Fig. |4fa) to the potential in 
Fig. I of the main text, we know that A s m A e and A s * A. 
Thus, we arrive to the results shown in the main text. 



